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Entanglement is often the key element in quantum information protocols. Here, we present 
schemes to generate robust photon entanglement in optomechanical interfaces via quantum interfer- 
ence. The schemes explore the excitation of the Bogoliubov dark mode and the destructive quantum 
interference between the bright modes of the interface, similar to electromagnetically induced trans- 
parency, to eliminate leading-order effects of the mechanical noise. Both continuous-variable and 
discrete-state entanglements that are robust against the mechanical noise can be achieved. 
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Introduction. Optomechanical systems can serve as key 
elements in hybrid quantum networks that connect opti- 
cal and microwave photons |l|, |2j . The mechanical modes 
can couple with cavity photons of distinct frequencies, as 
was demonstrated in recent experiments of optomechani- 
cally induced transparencies, normal mode splitting, cav- 
ity cooling approaching the ground state, and etc. [3-16]. 
Schemes for quantum state manipulation and transfer 
in optomechanical interfaces have been intensively stud- 
ied and the mechanical dark mode has recently been ob- 
served [30]. 

Entanglement is at the heart of many quantum infor- 
mation protocols. Various schemes for generating entan- 
glement in optomechanical systems, either between the 
photon states or between the cavity and the mechanical 
modes, have been studied [28 44], The amount of en- 
tanglement generated in these schemes is often limited 
by factors such as the stability of the system and the 
thermal noise of the mechanical mode. In stationary- 
state schemes where the cavities are continuously driven 
28L 2^ . 31-13], the stability conditions place an upper 
bound for the effective optomechanical coupling and con- 
strain the entanglement. In transient schemes where cav- 
ities are driven by pulses, the intrinsic nonlinearity of the 
radiation pressure force and th e thermal noise can con- 
strain the entanglement 
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Continuous- variable entanglement between cavity modes 
can be generated by the parametric amplifier Hamilto- 
nian H s = —g s {a\a2 + a\a\) with coupling g s , where a* 
(i = 1, 2) is the annihilation operator for mode i [45] . The 
cavity operators at time t can be written as cij(i) = f3%(r) 
in terms of the Bogoliubov modes 

/3i(r) = cosh(r)ai + ismh(r)a\ (la) 
&2(t) = cosh(r)a2 + isinh(r)a{ (lb) 

and the squeezing parameter r = g s t. When applied 
to the vacuum state IO1O2), this operation generates a 
two-mode squeezed vacuum state with the entanglement 
En = 2r log 2 (e), where the entanglement is quantified by 
logarithmic negativity [3]. In optomechanical systems, 
cavity modes couple with mechanical modes but do not 



couple directly with each other. The optomechanical cou- 
pling can induce mixing between the cavity and the me- 
chanical modes and expose the cavity modes to mechan- 
ical noise. In this work, we present schemes to generate 
strong photon entanglement that is robust against the 
mechanical thermal noise in an optomechanical interface 
by designing cavity operators in the time and frequency 
domains to have the Bogoliubov-like forms defined in 
Eqs. ljlallbp . In this interface, a Bogoliubov dark mode 
in the form of 0% (r) exists that does not contain mechan- 
ical operators. By combining the excitation of this mode 
and the destructive quantum interference between the 
bright modes of the interface, similar to electromagneti- 
cally induced transparency (EIT), the leading-order me- 
chanical components can be eliminated from the cavity 
operators. Robust entanglement can hence be achieved 
between cavity photons in the time domain and between 
cavity outputs in the frequency domain. Meanwhile, we 
show that robust entanglement in discrete photon states 
\ipen) = (|0i I2) ± |li02))/v / 2 can also be achieved in 
this interface. Our results demonstrate that optome- 
chanical systems can act as noise-resilient hubs in hybrid 
quantum networks to perform key quantum operations 
such as quantum state transfer and entanglement gen- 
eration. This facilitates the implementation of scalable 
hybrid quantum systems. Furthermore, the approaches 
used here can be extended to similar systems such as two 
cavity modes coupling with one noisy qubit to implement 
quantum operations. 

Quantum interface with Bogoliubov dark mode. The 
interface is composed of two cavity modes coupling 
with the same mechanical mode with the interaction 
^KGia\ai(b m + fr^J, where b m is the annihilation op- 
erator of the mechanical mode. One cavity is driven by 
red-detuned source with cavity detuning Ai to generate 
anti-Stokes processes and the other cavity is driven by 
blue-detuned source with cavity detuning A 2 to gener- 
ate Stokes processes as is illustrated in Fig|T](a,b). This 
model can be realized in many systems such as the hy- 
brid system of a microwave cavity and an optical cavity 
both coupling with a membrane 31- 34|. For general- 
ity, we use arbitrary units for the system parameters but 
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Figure 1: (a) The interface with three modes, (b) Cavity 
resonances at uj c i and driving sources at Udi with Ai = uidi — 
u) C i. (c) 5Ai (solid) and 8X2 (dashed) versus r for {k\,k<i) = 
(0.3, 0.2) (thick) and (0.2, 0.3) (thin) with g = 3 and j m = 
0.001 in arbitrary units. 



choose realistic relative values between these parameters. 
Let the detunings be — Ai = A2 = u> m with uj m being 
the mechanical frequency. After standard linearization, 
the effective Hamiltonian in the interaction picture of 



hujmb^bm can be written as 



Hi = %i(a}6 m + 6^,01) + ihg 2 (a\b\ n - a 2 b m ) (2) 



where g^s (i — 1,2) are the effective optomechanical cou- 
plings [47| . The environmental fluctuations can be repre- 
sented by the input operators a^* (i) for the cavities and 
binif) for the mechanical mode. The correlation functions 
for the input operators are (af^ (t)af^ (t 1 )) = 5(t—t') and 
{b ln {t)b\ n (t')) = {nth + l)S(t - t') in the high tempera- 
ture limit with thermal phonon number n t h [29]. The 
Langevin equation for this system is 



(3) 



idv(t)/dt = Mv(t) + WKv in (t) 



with v(t) — [cii(t), b m (t), al(t)] T for the system opera- 
tors, Vin(t) = [a$(t), b in (t), af^\t)] T for the input oper- 
ators, diagonal matrix K — Diag[/ci, 7 m , K2], and 



1 2 


9i 


) 




91 


1 2 


m 


(4) 





ig2 


-if J 





M 



where KiS (7m) are the cavity (mechanical) damping 
rates. With uj m S> <?i, Ki,^1%nth, the rotating wave ap- 
proximation has been applied in the above equations. 

The blue-detuned drive can induce instability in the in- 
terface, which affects the generation of entanglement. Us- 
ing the Routh-Hurwitz criterion [48j , we derive the stabil- 
ity conditions for this system which can be approximated 
as g\lg\ > max{K2/Ki, ki/k 2 } in the strong coupling 
regime with gi i2 3> Ki,7 m . This requires g\ > g 2 for the 
system to be stable. Thus, we can write g\ = go cosh(r) 
and g 2 = 9o sinh(r) with r = tanh -1 (52/91 )• 
Let at be an eigenmode of the matrix M with eigenvalue 
Ai. At zero damping, the eigenmodes are 



«i 



2 ; "2,3 



with eigenvalues Ai = and A2.3 = ±go- The mode ct\ 
only contains the cavity modes and we call it the Bogoli- 
ubov dark mode. The modes 02,3 contain both cavity and 
mechanical modes and we call them the bright modes. An 
interesting feature is the symmetry of the bring modes, 
which gives (a 2 + a^)/\/2 = /3i. At finite damping rates 
but in the strong coupling regime, the eigenmodes can be 
derived by treating the damping terms in the matrix M 
as perturbation. The eigenvectors of these modes form 
the matrix [53, E3] 



U 



-i sinh(r) 

Xi 

cosh(r) 



sh(r) 



X 2 



cosh(r) 

~7T~ 



x 2 



i sinh(r) 
V2 



(6) 



X4 



i sinh(r) 
V2 
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where xi = 0(Ki/go,jm/9a) is a first-order correction to 
the eigenmodes by the perturbation. The normalization 
conditions for the eigenvectors give UuUij = 5ij to 
the first order of the perturbation. We also have 

ai = 02 + xib m ; (0C2 + a 3 )/V2 = - V2x 3 b m (7) 

which contain first-order terms 0(xi)b m of the mechani- 
cal mode. The eigenvalues are also modified by first-order 
corrections with Ai = i5X\ and A 2 .3 = ±.g + i5X 2 . The 
small imaginary parts SXi are directly related to the sta- 
bility conditions of the interface which becomes unstable 
when one of SXi becomes positive. 

Robust entanglement in time domain. Under the adia- 
batic condition for the couplings [U 149] . the Langevin 
equation can be written in terms of the eigenmodes as 

ida{t)/dt = A(i)a(i) + iU^ 1 (t)-/Kv in (t) (8) 

where a — [a\ 1 a 2 , 0:3] T and A(t) — Diag [Ai, A2, A3] with 
a(t) = U T (t)v(t) and U T MU = A. The operators a(t) 
can be derived by integrating Eq.([5]). At zero damp- 
ing, ai(t) = ai(0) and a 2 ^(t) = exp(=Fz</j(i))a 2 ,3(0) with 
tp(t) = f Q dt'go(t'). Applying Eq.([S]), we derive that 



{t) = (0) cos <p(t) - ib m (0) sin <p(t) 



(9) 



which mixes the cavity and the mechanical modes, and 
P 2 {t) = ,$2(0). At time t n with ip(t n ) — mr for integer n, 
0(t n ) = (— l) n /3i(0). Hence, the Bogoliubov operators 
at time t n only contain the cavity modes ai(0) with the 
mechanical mode eliminated by the destructive quantum 
interference between the phase factors in 0:2,3 (t) [Hoj ] . 
For time-independent couplings with go and r, ip(t) = 
got. At t n = mr/go for even n, Eq.® gives 0{t n ) — 
0(0). Hence, ai(t n ) = ai(0) and the cavities return to 
their initial state. For odd n, 0(t n ) = — 0(Q). The 
cavity operators at t n can be derived as 



(3) 



ai(t n ) 
a\(t n ) 



cosh(2r) 
i sinh(2r) 



-i sinh(2r) 
cosh(2r) 



-ai(0) 

4(o) 



(10) 
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which gives a two-mode squeezed vacuum state with 
squeezing parameter 2r and entanglement 4r log 2 e when 
applied to the vacuum state. These operators do not con- 
tain the mechanical mode b m (0) and the entanglement at 
t n is thus not subject to the effect of thermal fluctuations 
in the initial mechanical state. We can also consider 
time-dependent couplings. For g\{t) — go cosh(Ai) and 
(t) = go sinh(Ai) under the adiabatic condition A <C go 
14 71 1 . we derive that at t n — mr/go for integer n, 



ttl(tn) 
4(*n) 



cosh(r) 
i sinh(r) 



-i sinh(r) 
cosh(r) 



(- 



l)"ai(0) 
4(0) 



(11) 

which describes a two-mode squeezed vacuum state with 
squeezing parameter r = Xt n and does not contain the 
mechanical component 6 m (0). 

At finite damping rates, both the thermal fluctuations in 
the initial mechanical state with average phonon number 
no and the input noise of the mechanical bath with ther- 
mal excitation number nth affect the entanglement gener- 
ated in the time-domain schemes. By solving Eq.©, the 
cavity operators di(t n ) can be derived [47], which include 
terms of 0(Ki/go)b m (0) and terms of 0(f dt 1 \f"fmbi n {t')) 
to leading order of m /go and Jm/go- We first assume 
no = nth- At t n — mr/go, the above terms affect the 
covariance matrix of the cavity state as 0(k? / 'gfynth 
and 0("f m nth/ go) respectively. While at any other time 
t ^ t n , the cavity operators contain the mechanical mode 
as O(l)b m (0) which affects the covariance matrix as nth- 
The destructive quantum interference at t n suppresses 
the thermal effects significantly by eliminating the lead- 
ing order terms of 6 m (0) from the cavity operators. The 
entanglement at t n is hence robust against the thermal 
noise. Note that other first-order terms in di(t) include 
0(ni/ go)(ii{0) due to the decay of the eigenmodes and 



0(J dt' 



Kid 



(i) 



(t')) due to cavity input noise, both of 



which affect the covariance matrix as O(ni/go)- 
The entanglement of the time-domain schemes is plotted 
in Fig[^a,b) [47j]. At finite temperature, sharp resonance 
peaks appear at t n , where the peak values decrease slowly 
with n t h but the peak widths decrease quickly with nth- 
This is because the cavity operators di(t n ± dt) at the 
small deviate St from t n contain the mechanical mode as 
O(go6t)b m (0) as can be derived from Eq.®. This term 
affects the covariance matrix and the entanglement as 
(goSt) nth ancl narrows the resonance peaks for large nth- 
In Figf^c), it is shown that the entanglement at the res- 
onant peaks remains sizable even for n t h ~ 10 4 , in sharp 
contrast to the stationary-state entanglement which de- 
creases to zero quickly. For the above, no is equal to nth- 
In Figf2fd), we study the effect of no on the entangle- 
ment. For no = 0, the entanglement increases with time. 
For finite no, the resonance peaks appear again and the 
peak widths decrease quickly with no- This result clearly 
verifies our analysis that the thermal fluctuations in ini- 
tial mechanical state narrow the resonances. 
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Figure 2: (Color online) Entanglement between cavity pho- 
tons, (a) and (b) En versus time for time-independent scheme 
[r — 1) and adiabatic scheme with r(t,2) = 1; no = nth = 
0, 10, 10 2 ,10 3 from top to bottom, (c) En versus nth for 
time- independent scheme at ti (solid), adiabatic scheme at t% 
(dashed), and stationary state (dash-dotted), (d) En for adi- 
abatic scheme at t% n t h = 10 3 and no = 0, 10, 10 2 , 10 3 from 
top to bottom, and no = n t h = (thin), (ki,kz) = (0.3,02) 
and other parameters are same as in FigfJ] 



Robust entanglement in frequency domain. For time- 
independent couplings, the Langevin equation can be 
written in the frequency domain as v(ui) — (Iuj — 
M)~ 1 \fKvi n (uj) with the frequency components defined 
byO(w) = f dtO(t)e iult /V27r. Using the input-output re- 
lation, the output operators can be derived as v out (uj) — 
S(u)v in {uS) with S(u) = (I - iVK(Iuj - M)- 1 ^). To 
study the entanglement between cavity outputs, we de- 
fine the operators dx\oj n ) — J dujg^(uj — Lu n )dx\uj) for 
Lj n = nAu (n integer) and x = in, out with the en- 
velop function £a(w) = 1/VAlj for u E 4f ) and 
3a(w) = otherwise. The commutation relations for 



these operators are [d x (oj m ),d x (u> n )] 



Smn^ij, which 



ensures that the covariance matrix of their quadrature 
variables can be directly used to calculate entanglement 
(iH , 47 1 . In Figl2Ia,b), entanglement between the out- 
put modes is plotted versus u n . Three resonance peaks 
appear at uj n = 0, ±<?o, corresponding to the three eigen- 
modes cti respectively. At ljo [i.e. u) n for n = 0), 
the entanglement decreases slowly at large phonon num- 
ber nth an d is robust against the thermal noise even 
for nth = 10 4 . While at uj n = ±<7oj the entanglement 
decreases quickly to reach zero. This is illustrated in 
FigHcd). 

For input noise at frequency uj n , the excitations of eigen- 
modes are a(cj n ) = i(Iw n — A)~ 1 £/ T \/~Kvi n (uj n ) [47]]. At 
wo, the Bogoliubov dark mode is strongly excited as 

• cosh(r) \ 



ai(wo) = ( 



sinh(r) 
<5Ai 



J <5Ai 



<5Ai 



Kv in (w Q ) (12) 



with x\/5\\ ~ O(l/<7o)- The bright modes are weakly 
excited as «2,3(wo) oc \/go with the relation 



(a 2 (wo) + a 3 {u>o)) / V2 - 



(vo)/go- (13) 
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Figure 3: (Color online) Entanglement between output pho- 
tons, (a) and (b) En versus to n for («i,/C2) = (0.3,0.2) and 
(0.2,0.3); n th = 0, 10, 10 2 , 10 3 from top to bottom, (c) and 
(d) E N versus n th for (ki,k 2 ) = (0.3,0.2) and (0.2,0.3) at 
cjo (solid), u) n = ±<?o (dashed), and stationary state (dash- 
dotted). Other parameters are same as in Fig[T] 



The cavity operators a,i{u>o) can be derived from these 
results and contain the cavity inputs as Q(a^ (^n) / 
and the mechanical input as 0(^/7 m &m(wo) / go) |47| . The 
mechanical input is a factor Ki /go smaller than the cavity 
inputs when 'yrn.nth ~ This is due to the destructive 
quantum interference between the excitations of 0:2(^0) 
and a^(u>o)- As a result, the thermal effect in the co- 
variance matrix of the output operators a^ t (u>o) is sup- 
pressed by a factor of {ni/go) 2 and the entanglement at 
wo is robust against the mechanical thermal noise. At 
= go, ^1,3(90) ~ I/50 weakly excited. While one 
of the bright modes 0:2(30) ~ 1/6X2 is strongly excited 
due to its resonance with the input noise and contains 
the mechanical input as O(y/j m bi n (go) / 6X2). This large 

mechanical term cannot be eliminated from a^ t (go) an d 
will impair the entanglement quickly as n t h increases. 
Similar results can be derived for u n = —go- 

The entanglement here depends strongly on the cavity 
damping rates as is shown in Figj3ja,b). This is due 
to the dependence of 6Xi on the damping rates. In the 
stable regime, 6X\ } 2 < 0. For K4 > k 2 , as the squeezing 
parameter r increases towards the unstable regime, \6Xi\ 
becomes much smaller than | <5 A2 1 with <5Ai — > as is plot- 
ted in FigQJc). In this regime, the excitation of the Bo- 
goliubov dark mode at loo becomes significantly stronger 
than the excitation of the bright modes at u n — ±go so 
that the entanglement at coo is stronger. For K2 > K\, 
as r increases towards the unstable regime, 6X2 — > 
and the entanglement at uj n — ±go is stronger. Mean- 
while, regardless of the damping rates, the entanglement 
at ujq is always robust against the thermal noise. This 
tells us that in a hybrid interface with very different cav- 
ity damping rates, by applying red-detuned drive to the 
cavity mode with larger damping rate and applying blue- 



detuned drive to the other cavity to get n\ > K2, stronger 
entanglement that is robust against thermal noise can be 
achieved. 

Robust entanglement in discrete states. This interface 
can also be used to generate entanglement in discrete 
states such as (|li0 2 ) ± |0il 2 >) /a/2 [51] . Let the two 
cavities both be driven by red-detuned sources with 
— Aj = ui m . In (22I [23|], this setup was studied for high 
fidelity transfer of quantum states. Assume the couplings 
to be gi(t) = gosin(At) and 02 (i) = — Socos(At), varying 
adiabatically with X < go- At if = 7r/4A and by choosing 
A = go/4n for integer n, the cavity operators become 

/ ai(t f ) \ _ 1_ ( 1 -1\( oi(0) \ ( 

Ww~ V2 V 1 1 A(- 1 ) n °2(°)y 

It can be proven that for an initial cavity state | I1O2) > the 
final state of the cavities is \ip en ) = (II1O2) + |0il2))/v / 2- 
Similarly, for the initial state 1 0i I2) , the final state is 
\ipen) = (| li0 2 ) - |0il 2 ))/v / 2. The effect of the thermal 
noise can be studied by solving the Langevin equation us- 
ing perturbation theory [47]. The cavity operators aj(i/) 
contain the mechanical mode as 0(Kj/ go)b m (0) which is 
suppressed by factor Hi /go due to the destructive inter- 
ference. Hence, the discrete-state entanglement is also 
robust against thermal noise. 

Conclusions. We study an optomechanical interface for 
the generation of photon entanglement that is robust 
against the mechanical noise in both the time and the 
frequency domains. Due to the excitation of the Bogoli- 
ubov dark mode and the quantum interference between 
the bright modes, the effect of the mechanical noise is 
significantly suppressed. Both continuous-variable and 
discrete-state entanglements can be generated. When 
combined with the state transfer schemes, this quantum 
interface provides a promising building block for hybrid 
quantum networks and for quantum state engineering. 
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